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The geometrically-frustrated, triangular antiferromagnet GGG exhibits a rich mix of short-range 
order and isolated quantum states. We investigate the effects of up to 1% Neodymium substitution 
for Gallium on the ac magnetic response at temperatures below 1 K in both the linear and nonlinear 
regimes. Substitutional disorder actually drives the system towards a more perfectly frustrated 
state, apparently compensating for the effect of imperfect Gadolinium/Gallium stoichiometry, while 
at the same time more closely demarcating the boundaries of isolated, coherent clusters composed 
of hundreds of spins. Optical measurements of the local Nd environment substantiate the picture 
of an increased frustration index with doping. 

PACS numbers: 75.50.Lk, 75.40.Gb, 75.45.+j, 75.30.Hx 


Geometrically frustrated magnets suppress the ability 
of spins to freeze to temperatures well below their nat¬ 
ural ordering scales and, in so doing, open to experi¬ 
mental scrutiny low temperature fluctuations, extended 
quantum states, and thermodynamically isolated degrees 
of freedom. The topology of the interspin interactions 
makes it impossible to simultaneously minimize all of 
the pairwise interaction energies^ and, in principle, a 
fully frustrated system cannot form a long-range ordered 
state, even at absolute zero. Short-range effects can dom¬ 
inate the macroscopic behavior, with self-organization of 
spins into such quantities as quantum protectorates: co¬ 
herent states that are topologically isolated from many 
sources of environmental decoherence PS The introduc¬ 
tion of structural or compositional disorder to such a 
frustrated system has the potential to tune the frus¬ 
tration effects and their macroscopic expressions. It is 
not clear a priori , however, whether additional disor¬ 
der enhances self-organization by further isolating the 
spin clusters, or actually relieves frustration by break¬ 
ing the symmetry of the interspin couplings and hence 
degrading the isolation bet ween any extended quantum 
states and the environment!^ Here, we study the ef¬ 
fects of progressively increasing the degree of disorder in 
a geometrically frustrated antiferromagnet by introduc¬ 
ing controlled amounts of a magnetic dopant. 

The parent compound, Gadolinium Gallium Garnet 
(Gd 3 GasOi 2 or GGG), has a cubic lattice in which the 
magnetic Gd 3+ ions are located on two interpenetrating 
sub-lattices composed of corner sharing triangles with a 
Heisenberg spin symmetry and a single ion anisotropy 
of less than 0.04 kP Nearest neighbors are connected 
by corner sharing triangles on each sub-lattice, coupled 
by a 1.5 K antiferromagnetic exchange interaction and a 
0.7 K dipolar interaction, and the two sub-lattices inter¬ 
act via next-nearest-neighbor interactions. The resulting 
Weiss temperature is of order 2.3 K. 7 The combination 
of the antiferromagnetic interaction and the topology 
of couplings provides a macroscopic number of equally- 


probable spin configurations at absolute zero. 8 

Several experimental studies of GGG have probed the 
subtle nature of the ground state, with differing re¬ 
sults. Early bulk thermodynamic measurements 9 found 
evidence of spin freezing and an associated spin glass 
transition at T ~ 130 mK. By contrast, muon spin 
relaxatioifhl] 1: 0 and Mossbauei^ measurements found 
that the spins remained unfrozen to temperatures as low 
as 25 mK. Powder neutron measurements observed short- 
range correlations at elevated temperatures followed by 
a 140 mK transition to a state with longer range cor¬ 
relations and spin-liquid type behavior. 13 Subsequently, 
a combination of linear and nonlinear ac susceptibility 
measurements found a state consisting of a set of “quan¬ 
tum protectorates”, extended coherent states that are 
topologically decoupled from the background, coexisting 
with at least short-range antiferromagnetic order onset- 
ting at T ~ 100 rnKpdl More recent high-resolution heat 
capacity measurements 15 did not observe any evidence 
of a spin glass transition in the 120-200 mK range. The 
interesting spin dynamics all occur at temperatures be¬ 
low 200 mK, over an order of magnitude lower than the 
Weiss temperature, confirming the picture of a highly 
frustrated spin configuration. 

One effect contributing to this range of results is the 
difficulty in growing purely stoichiometric crystals of 
GGG, as excess magnetic Gd 3+ ions tend to randomly 
substitute for nonmagnetic Ga 3+ ions! 16 ! 17 ] Gadolinium 
ions occupying gallium sites have a nearest-magnetic- 
neighbor distance that is 9% smaller than the normal Gd- 
Gd distance in GGG, resulting in both lattice strain and 
a local increase in the magnetic interaction strength. 15 
Against this backdrop, we examine a series of GGG crys¬ 
tals with up to 1% of Nd 3+ dopants substituting for Gd 3+ 
ions. The Nd ions introduce a controlled degree of disor¬ 
der into the magnetic lattice that breaks the local symme¬ 
try of the magnetic interactions. By controlling the level 
of this quenched disorder, we can study the role it plays 
in the formation and evolution of the low-temperature 
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FIG. 1. Imaginary part of the ac magnetic susceptibility x" 
as a function of frequency / at a series of temperatures T 
for GGGiNd^ with x = 0, 0.1%, and 1%. The susceptibility 
was measured using a 40 nrOe probe field to ensure linear re¬ 
sponse. A low-frequency mode at 50 Hz and a high-frequency 
mode that peaks above the 10 kHz measurement window char¬ 
acterize the global response of all three materials. For x = 0 
and 0.1%, a flat low-frequency response develops at 85 to 95 
mK, corresponding to fluctuations associated with ordering. 
At x = 1.0%, this feature is not observed, indicating that any 
ordering is suppressed below 30 mK. 


magnetic phases. 

We performed ac susceptibility measurements on a se¬ 
ries of (5 x 5 x 10) mm 3 single crystals of GGG:Nd a ., 
with x = 0, 0.1, and 1.0% (Princeton Scientific). X- 
ray lattice constant measurements showed that all three 
crystals exhibited a similar degree of off-stoichiometry of 
approximately 3%. We employed a gradiometric suscep- 
tometer attached to the cold finger of a helium dilution 
refrigerator, using sapphire rods pressed against the faces 
of the crystals to provide a thermal link to the cryostat. 
The complex ac susceptibility was measured in the linear- 
response limit using 40 mOe of probe field and mapped as 
functions of frequency (1 Hz to 10 kHz) and temperature 
(30 to 300 mK) for all three dopings. 

The nonlinear response of the system provides the most 
acute insight into the behavior of the local, isolated spin 
clusters. To that end, we measured the susceptibility as 
a function of drive amplitude (mOe to Oe) at fixed fre¬ 
quency. In addition to the direct effect on the GGG 
crystals, large ac magnetic fields, particularly at kilo¬ 
hertz frequencies, also induce eddy currents in metallic 
mounts, resulting in significant Ohmic heating. To al¬ 
leviate this issue, we mounted the susceptometer on a 
carbon-fiber/sapphire framework 18 designed to minimize 
eddy currents while maintaining a high thermal conduc¬ 
tivity for efficient thermal linkage to the cryostat. The 
ac probe field was applied using a duty cycle of 1-2% to 
further reduce the heat load. 
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FIG. 2. Real part of the susceptibility as a function of 
temperature T at a series of frequencies / for GGG:Nda, with 
x = 0, 0.1%, and 1%. The functional form and tempera¬ 
ture scale for all three doping levels is similar; however, the 
static incoherent background is suppressed for x = 1.0% as 
frustration is enhanced by doping. 

Finally, we exploited the coupling of the Nd ions to 
the Gd ions to probe optically the local magnetic environ¬ 
ment. Lifetime data was calculated from the homogenous 
linewidth of the Ri —► Z 3 transition of the 4 F 3 / 2 —► 4 /g/ 2 
ground to first excited state emission multiplets of the 
Nd 3+ ion. The samples were mounted in a liquid helium 
flow cryostat and excited with a Ti:Sapphire laser tuned 
to 808 nm. The excitation wavelength was chosen to be 
as close to the absorption band edge as possible for the 
4 -^ 3/2 —t 4 dg /2 emission lines. The emission is collected 
in reflection and dispersed by an Acton 300i spectrome¬ 
ter onto a tlrermoelectrically cooled CCD with spectral 
resolution of 0.18 nm. 

The complex ac susceptibility at audio frequencies de¬ 
lineates the nature and timescales of the low-energy mag¬ 
netic modes of the material. Peaks in the imaginary 
component of the susceptibility, x" 1 yield information on 
characteristic relaxation times. The presence of a flat re¬ 
sponse in the / —> 0 limit corresponds to 1 // noise in the 
magnetization, a signature of a spin state with sca le in- 
variance and hence a complex hierarchy of energies! 19 ! 20 ! 
We plot in Fig. [T] imaginary susceptibility spectra for 
GGGrNdz (x = 0, 0.1, 1%) for 30 < T < 300 mK. For all 
three samples, we observe two peaks in the susceptibility, 
one at 50 Hz and one that appears to lie above the 10 
kHz upper limit of the measurement. This corresponds 
to two characteristic relaxation timescales in the system, 
one at 20 msec and one shorter than 100 /xsec. 

The characteristic long-time relaxation at 50 Hz does 
not depend strongly on either T or x. By contrast, the 
behavior of the system in the low-frequency limit changes 
markedly with x. For pure GGG and GGG with 0.1% 
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Nd, we observe a clear / —> 0 plateau in the imaginary 
susceptibility at T ~ 85 to 95 mK. No such plateau ap¬ 
pears in the 1.0% doped material. The appearance of 
such a plateau is often correlated with the entry into 
a spin glass stated 21 If that were the case here, the 
plateau should remain robust for all T below the glass 
transition temperature. Instead, the plateau behavior is 
only stable over a finite temperature band, below which 
the susceptibility approaches zero with finite slope. We 
associate the plateau with antiferromagnetic ordering, at 
least on the 100 A scaled The introduction of Nd sup¬ 
presses the ordering to T < 30 mK, effectively restoring 
a higher degree of geometric frustration. 

The relative frustration index as a function of doping 
also can be discerned by looking at the real component, 
x', of the linear susceptibility, as shown in Fig. [2] for all 
three sample concentrations, x' consists of a frequency- 
dependent peaked component riding on a frequency inde¬ 
pendent background. The peak reflects the freezing out 
of magnetic modes at the pertinent time scale of the in¬ 
terrogation frequency, but conspicuously does not go to 
zero (note the suppressed zero of Fig. [2]). We attribute 
this background to a Heisenberg spin bath made avail¬ 
able by imperfections in the crystal and deviations from 
stoichiometry. Rather than increasing the background 
susceptibility, the Nd ions at the highest level of substi¬ 
tution effectively reduce the degree of disorder and push 
the system towards a more perfectly frustrated state. 

We can ask whether this apparently counter-intuitive 
result also applies to the quantum protectorates previ¬ 
ously reported for undoped GGGP^ An experimental sig¬ 
nature of these coherent spin clusters in GGG occurs 
in the nonlinear susceptibility. As a function of drive 
field amplitude, the magnetization traces out a Brillouin 
function, with approach to saturation above 1 Oe. This 
behavior is characteristic of clusters that act as large ef¬ 
fective spins excited between quantized statesPf Similar 
field-induced formation of extended states has been ob¬ 
served as well in the dilute Ising magnet LiHo a; Yi_ a: F 4 , 
where quantum fluctuations rather than geometric frus¬ 
tration preserve the isolated degrees of freedomP^HMI 
Here we explore the effects of the Nd doping on these 
cluster states in GGG. 

We plot in Fig. [3] the nonlinear susceptibility at / = 5 
kHz for the three dopant concentrations over almost a 
decade in T. The response at / = 15 Hz is qualita¬ 
tively similar. There are two distinct regimes in tem¬ 
perature. For T > 175 mK, the susceptibility is essen¬ 
tially constant at low applied field, followed by a con¬ 
tinuous drop once the applied field is above a threshold 
value. The shape of the susceptibility above the thresh¬ 
old can be described by the derivative of the Brillouin 
function, x(if) = N c N s HBgJ + const., where 
V = Ns k B B T H ■ B J = ^coth (HAly) - ajcoth (±y), 
and for GGG J = 7/2 and g = 2. N s is the number of mi¬ 
croscopic spins bound in an individual cluster, N c is the 
total number of clusters per unit volume, and the con¬ 
stant term accounts for the susceptibility of spins that 
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FIG. 3. Real part of the susceptibility x! of GGGiNd^ in 
the nonlinear response regime, measured at / = 5 kHz as a 
function of the applied probe field H at a series of tempera¬ 
tures. At elevated T, there is clear threshold behavior in the 
susceptibility and the nonlinear response can be fit to a Bril¬ 
louin function characteristic of coherent clusters with a large 
effective spin (see text for details). Below T = 175 mK, the 
cluster dynamics freeze out and the small increase in \ arises 
from inductive heating of the susceptometer (and hence the 
sample) at large H. 


are not bound up in clusters. Even in the presence of 
doping-induced disorder, clusters acting as single large 
effective spins continue to form. For T < 160 mK, the 
spin clusters freeze out and the nonlinear contribution to 
the susceptibility disappears. The slow rise in x(U) at 
the highest drive fields is most likely due to Ohmic heat¬ 
ing that is at worst 25 mK for a 2.75 Oe/5 kHz applied 
field at lowest T by reference to the linear-response data 
of Fig. [2] 

The fit to the Brillouin form for T > 175 mK helps 
specify the nature of the cluster dynamics. We show in 
Fig. a the aggregate moment of the spins bound up in 
coherent clusters at / = 15 Hz, as well as the threshold 
activation field for the onset of the nonlinear response as 
a function of T. The typical cluster size is 150-200 spins, 
corresponding to 6 to 8 unit cells of the GGG lattice. The 
moment of the clusters only depends weakly on T, until 
the entire mode freezes out (Fig. [3]). There is, however, 
a definite compositional dependence (Fig. The 0 

and 0.1% Nd-doped samples have comparable moments, 
whereas the 1% crystal has a net cluster moment reduced 
by 20%. The Nd does not act to improve the cluster cor¬ 
relations (by, e.g., acting as nucleation sites for coherent 
clusters), but instead serves to reduce the disorder in the 
GGG lattice and more closely demarcates the boundaries 
of the isolated spin clusters. 

While the net moment of the clusters does not vary 
appreciably with T, the threshold field required to excite 
the coherent response is thermally activated and obeys an 
Arrhenius form (Fig. |4Jd) with an energy barrier height 
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that spin-spin interactions are the dominant effect. As 



FIG. 4. Cluster dynamics of GGGiNd^ at / = 15 Hz. (a) 
Aggregate moment of all spins involved in coherent clusters 
as a function of temperature T for x = 0, 0.1, and 1.0%. 
Dashed lines show the mean values for the three doping con¬ 
centrations. Increased doping reduces the cluster size, (b) 
Activation field vs. inverse temperature. The Arrhenius form 
establishes the germane energy scales for cluster spin flips as 
a function of x. 


of 0.48 ± 0.05 K. The energy gJ^sH to flip a cluster of 
200 spins with J = 7/2, g = 2 and a threshold field for 
the nonlinear response ~ 1 Oe is approximately 0.2 K, 
consonant with the measured barrier height and the scale 
set by the dipolar coupling. Although the barrier height 
is equal within error bars for all three Nd concentrations, 
the threshold field for the 1% doped sample is reduced 
by 25% in comparison to the more lightly doped samples 
as the smaller clusters more easily change orientation in 
response to the driving field. 

We show in Fig. [5] the decay time of the optical Ri — ► 
Z$ emission line of Nd 3+ in GGG:Nd x as a function of 
temperature for x = 0.1 and 1%. Typically, as a system 
is cooled, the homogeneous linewidth (inversely propor¬ 
tional to the decay time) of a dipole transition decreases 
with the suppression of the phonon density of states. The 
contrary trend observed here indicates the availability of 
a different decoherence mechanism, and while this de¬ 
creasing trend in decohere nce is not seen in other Nd 3+ 
garnets like YAG and YAPp^H p has been reported in 
Mn 4+ doped GGG, where it was attributed to spin-spin 
coupling between the Gd and the dopant ions! 2 ' It is im¬ 
portant to note that that the excited state lifetime is long 
in rare earth doped crystals.^ and thus, judging by the 
extremely fast decay time seen in the data, we can infer 


FIG. 5. Decay time of the R\ — > Z$ emission line of Nd 3+ 
in GGG:Nda; as a function of temperature for x = 0.1% and 
1.0%. The faster decoherence time with increased dopant con¬ 
centration supports the notion of increased frustration.Lines 
are guides to the eye. 

the Nd 3+ ions fill in the Gd vacancies in the GGG lat¬ 
tice, the local environment becomes less disordered, the 
frustration index increases, and the decoherence time de¬ 
creases. 

Using a Nd doping series to vary the effective degree 
of disorder in GGG, we demonstrate that increasing the 
Nd concentration actually reduces the effects of disorder 
and enhances the signatures of frustration. We posit that 
the Nd compensates for the disorder arising from the off 
stoichiometric growth process and, in so doing, permits 
the intrinsic geometric frustration of the lattice to sup¬ 
press the formation of an ordered state. The nonlinear 
susceptibility reveals quantitative aspects of the coher¬ 
ent spin clusters and their tenability with disorder and 
frustration. In principle, the Nd ions could be marshaled 
to interrogate locally the isolated cluster dynamics via 
optical studies at lower T. 
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